Off-axis digital hologram reconstruction: some practical 

considerations 



N. Verrier and M. Atlan 

Institut Langevin. Fondation Pierre- Gilles de Gennes. Centre National de la Recherche Scientifique 
(CNRS) UMR 7587, Institut National de la Sante et de la Recherche Medicale (INSERM) U 979, 
Universite Pierre et Marie Curie (UPMC), Universite Paris 7. Ecole Superieure de Physique et de Chimie 
Industrielles (ESP CI ParisTech) -10 rue Vauquelin. 75005 Paris. France 

nicolas. verrier@espci.fr 
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some of the main numerical processing methods, based either on the Fourier transform 
interpretation of the propagation integral or on its linear system counterpart, is reported. 
Less common methods such as adjustable magnification reconstruction schemes or Fresnelet 
decomposition are presented and applied to the digital treatment of off- axis holograms. 
The influence of experimental parameters on the classical hologram reconstruction methods 
is assessed, offering guidelines for optimal image rendering regarding to the hologram 
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1. Introduction 

Optical holography consists of the acquisition of 
images from diffracted optical field measurements. 
Holographic imaging was initially proposed by 
Gabor [2j for electron microscopy. Holograms were 
recorded on high resolution photographic plates. Due 
to the subsequently-called "in-line" configuration, 
holograms were stained with twin-image and zero- 
order contributions, which overlapped with the signal 
image ^J. Originally recorded with the red radiation 
of Mercury lamps, holograms were increasingly 
recorded with laser light sources, which gave much 
more reliable results. In 1962, Leith and Upatnieks 
proposed to introduce an off-axis reference beam |3] 
to separate, in the spatial frequency domain, the 
real image from the twin-image and zero-order 
diffraction terms. However, holograms were still to 
be reconstructed by optical means. 

The first digital reconstructions of optically- 



measured holograms were realized by Goodman [3] 
and further by Kronrod [5j (in Russian. More 
details can be found in Ref. [6j) in the early 1970's. 
Here, optically magnified parts of the holograms 
are digitally sampled and then reconstructed using 
Fourier-transform based routines. Digit alizat ion of 
optical holograms allowed, for instance, to improve 
reconstruction quality [7J, to retrieve information 
about the phase of the recording wave [8l[9], and to 
treat holograms without reconstruction [10]. One of 
the major breakthroughs in holographic imaging was 
initiated, by Schnars, with direct recording of digital 
holograms [llj. Charge-Coupled Device (CCD) 
and Complementary Metal Oxide Semiconductor 
(CMOS) digital sensor arrays enabled the acqui- 
sition and numerical processing of high-resolution 
holograms at fast rates. 

Intrinsic properties of holographic imaging allow 
this technique to be used in a wide range of domains 
such as fluid mechanics [T2H2Q] . biomedical imag- 
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ing [2IH2S1 or mechanical vibration analysis P7H55] . 
Democratization of high resolution CCD and CMOS 
sensors played a major role in the development 
of digital reconstruction techniques. For instance, 
reconstruction of phase-only [36j, shifted [37J, tilted 
or aberrated data [3MiD] has been successfully 
demonstrated. Inverse-problem approaches make 
it possible to improve object localization and field 
of view in the reconstructed hologram [4T| |42]. 
Compressive sensing based approaches are also to 
be considered when working in noisy or low-light 
conditions [131 HI]. Moreover, owing to the massive 
parallelization of image processing calculations 
by Graphics Processing Units (GPU), hologram 
reconstruction can be performed in real-time [45H48] . 

In this paper, we will describe most of the common 
off- axis digital holographic reconstruction schemes, 
and discuss their applicability. After some brief re- 
minders about digital holographic recording, we will 
present the main reconstruction approaches, involv- 
ing one to three Fourier transforms. Then a discussion 
about reconstruction with adjustable magnification is 
proposed. Methods to tackle aliases and replicas are 
proposed, leading to high quality magnified recon- 
structions. The use of Fresnelet transform will also be 
discussed. Reconstruction methods will be assessed 
experimentally with optically-acquired off-axis holo- 
grams, to provide insight into their respective suit- 
ability towards targeted applications. 

2. Fresnel holography bases 

Digital holography typically consists of recording an 
optical field emerging from an illuminated object 
in a diffraction plane {e.g. in free-space propaga- 
tion conditions), and numerically calculating, from 
diffraction models, the field distribution in the re- 
construction plane. In practice, optical holograms are 
measured-out from the interference of the diffracted 
beam beating against a reference beam, which is 
not disturbed by the object to be analyzed. One 
of the object-reference cross terms typically yields 
a complex- valued map {i.e. quadrature-resolved : in 
amplitude and phase) of the diffraction field in the 
sensor plane. The complex- valued measurement con- 



tains relevant information about the local retarda- 
tion of the diffracted field. Phase-shifting [49] and 
frequency-shifting [50j techniques were proposed to 
record the diffraction field in quadrature. The inter- 
ference pattern, recorded by sensor array, can be ex- 
pressed as [5T] 

E{x,y) = \n{x,y)f^\0{x,y)\^ 

+ O* (x, y) n (x, y)^0 (x, y) 7^* (x, y) , (1) 

where IZ and O denote reference and object optical 
fields respectively. Starred (*) symbols are associated 
with complex conjugate values. 



Obj. arm 




(a) (b) 



Fig. 1: (a) Hologram recording in off-axis configura- 
tion, (b) Spatial frequency representation of off-axis 
holograms. 

The interference between reference and object 
beams can be recorded within a wide range of con- 
figurations. These can be grouped in two main cate- 
gories: in-line [51J and off-axis configurations [55]. For 
our applications, off-axis holograms will be recorded 
with a Mach-Zehnder configuration. The final part 
of our off- axis Mach-Zehnder configuration is illus- 
trated on Fig. ([T|) (a). Here, reference and object 
beam are combined, with a relative angle a, using 
a non-polarizing beam-splitting/combining cube. It 
should be noted that a should be chosen so as to ful- 
fill the sampling theorem. The maximal value, lead- 
ing to a correct sampling of the interference pattern 
is therefore given, under paraxial conditions, by: 

A 

where Ax denotes the sampling rate of the recording 
device. 
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This off-axis angle results in separation, of the four 
terms of Eq. ([T]), in the spatial frequency domain. 
This aspect is proposed Fig. ([1]) (b). The central part 
of the hologram spectrum {\7l {x^y)\^ -\- \ {x,y)f) is 
known as the autocorrelation term, its size is associ- 
ated with the highest spatial frequencies of the object, 
denoted by B. Real and twin images of the object are 
respectively given by OIZ^ and T^O* . These two terms 
are twice as small as the autocorrelation term. To im- 
prove the reconstruction quality, autocorrelation and 
twin image terms have to be canceled. This can be 
achieved either by spatial filtering [53J or phase shift- 
ing [54J. Thus doing makes it possible to reconstruct 
the real image term only. 

In the following part, we will focus on hologram re- 
construction. After a brief reminder about the Fresnel 
transform, we will discuss its main digital implemen- 
tations. 

3. Digital hologram reconstruction 

Digital reconstruction of an hologram consists in a a 
posteriori refocusing over the original object, which 
can be performed by calculating backward propaga- 
tion of the light from the hologram to the reconstruc- 
tion plane. This process is equivalent to positioning 
the recorded hologram back into the reference beam. 
Reference beam therefore become the reconstruction 
beam. Using the Huygens-Fresnel principle, one can 
infer an integral formulation of the intensity Ej-ec in 
the reconstruction plane, from an off axis recorded 
hologram E ^5] 

Erec (e, V) = -4 / E y) ^^^^^^dxdy. (3) 
A Jr2 r 

Here, the distance r is given by: 



and Eq. (|3|) is rewritten as: 



(4) 



(x, y) and ((f , 77) denote the spatial coordinates in the 
hologram and reconstruction plane respectively. It 
should be noted that, the hologram E (x, y) has been 
recorded an off-axis configuration. Under Fresnel ap- 
proximation, when >> ^ {£, 
r can be approximated by: 



xf ^{r]-yf 



1 2 



1 



y-T] 



(5) 



exp (- 



iXz 



E{x,y) 



^^[(^-0 ^{y-v) \jdxdy. (6) 

This relationship will be used in the remainder of this 
paper to perform reconstruction of off-axis intensity 
holograms. 

As far as variables in Eq. ([6]) are separable, all the 
discrete formulations will be derived in the ID case. 
Generalization in two dimensions is straightforward. 
The ID discrete Fresnel transform is defined by: 



Erec {P) 



exp [^Z 



N-1 



iXz 



exp (i^p'^Ae) ^^(n) 



n=0 



X exp (^i^n'^Ax'^^ exp ^-i^npAxA^^ , (7) 

where nAx and pA^ respectively denote the spatial 
coordinate in the CCD and reconstruction plane, and 
N is the number of sampling points. 

Direct implementation of Eq. ([7j) is a time consum- 
ing process. Starting from Eq. (|6]), one can realize 
that efficient computational schemes can be designed 
to implement digital holographic reconstruction. This 
makes it possible to separate reconstruction methods 
into two main families: the Fourier based approaches 
(based on the use of a single fast Fourier transform 
(FFT)) well suited for imaging extended objects lo- 
calized far from the CCD or CMOS sensor, and the 
convolution methods, computed by using two or three 
FFTs. These methods are well adapted for the recon- 
struction of holograms, of small lateral dimensions, 
recorded near the imaging device. Alternative meth- 
ods can be considered when an adjustable magnifica- 
tion or advanced filtering techniques are needed. 

In the remainder of this section, we will detail the 
different computational approaches and apply these 
to the reconstruction of digital holograms. 

3. A . Single- FFT method 

Efficient implementation of Eq. ([7]) can be performed 



using FFT algorithm | 



]. In this case pixel pitches 



in both reconstruction (A^) and CCD plane (Ax) are 
related by: 

A? = T^^. (8) 
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Therefore, Eq. ([7]) can be rewritten as: 



N-1 



exp(^z) ____ ^. nXzp^ \ 



■ exp z 



iXz ^^^^ V N'^^x'^ , 
E (n) exp (^i^n^Ax^^ exp (^-i27r^^ . (9) 

This relationship is therefore easily computed by 



^rec (0 



X J" ^yE {x) exp 



TT 



■ (10) 



It should be noted that, within this configuration, the 
ratio between the reconstructed horizon and the sen- 
sor array extension (which can be abusively denoted 
as the magnification of the reconstruction method) 
is closely linked to the reconstruction distance i.e. 
7 = A^/Ax = \z/{NAx^). In the remainder of this 
paper, the intrinsic magnification of the single-FFT 
implementation of the reconstruction integral will be 
denoted by 70 = \z/{NAx'^). 

3.B. Convolution based approaches 

Holographic reconstruction can be viewed as a linear 
system. As matter of fact, Eq. (j6]) is the mathemati- 
cal expression of the spatial convolution between the 
hologram, and the Fresnel impulse response function 
hz^ which is defined by (omitting the multiplicative 
constant): 

(11) 



hz {x) = exp {^~^^^^ 



Convolution based approaches lead to unitary mag- 
nification, namely = Ax. 

3.B.I. "Three-FFT algorithm" 

Computation of the convolution product between the 
hologram and the holographic impulse response can 
be efficiently implemented in Fourier domain. Using 
fast Fourier transform algorithms, Eq. (|6|) can be 
computed as: 



^rec (0 



exp {^z) 



iXz 



J^-'[T{E{x)}T{hz (x)}], 

(12) 

where T and respectively stand for the Fourier 
transform and its inverse. 



3.B.2. Angular spectrum propagation 

This method is based on the propagation of the angu- 
lar spectrum of the hologram. The angular spectrum 
transfer function is given by (58) : 

^.^z I _ 1 , 
2i- 



H {u) ^ exp 



A 



- \2 2 

-Xu 



(13) 



where u is the spatial frequency in Fourier domain. 
Using Eq. (p!3|) . hologram reconstruction can be per- 
formed: 



^rec (0 



^^F-'[F{E{x)}H{u)], (14) 
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Fig. 2: Angular acceptance of digital holographic re- 
construction process. Solid lines are associated with 
the 1-FFT reconstruction, and dashed lines corre- 
spond to the convolution approaches. 

S.C. Algorithms with adjustable magnification 

Neither the single-FFT approach nor the convolution 
based methods allow the adjustment the magnifica- 
tion of the reconstructed hologram. As a matter of 
fact, in a single-FFT hologram processing, magnifi- 
cation depends on the recording wavelength and dis- 
tance, whereas it remains constant using a convolu- 
tion approach. In the latter case, the magnification 
is unitary (A^ = Ax). This aspect is illustrated by 
fig. ([2]). Here, the evolution of the reconstructed hori- 
zon is represented with respect to the reconstruction 
distance. The solid lines are associated with the 1- 
FFT reconstruction scheme, and the doted lines are 
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the reconstruction horizon of the convolution based 
reconstruction approaches. It should be noted that 
for z = NAx'^/X (this distance is determined by tak- 
ing Ax = in Eq. (|8])), 1-FFT and convolution 
approaches exhibit the same magnification. 

Working with an adjustable magnification algo- 
rithm is a great opportunity to make the recon- 
struction horizon independent from hologram record- 
ing parameters. Domains such as multi-wavelength 
holography benefit from this property [591160]. 

Several approaches have been proposed to allow 
magnification adjustment. Ferraro used zero-padding 
to control the reconstructed horizon and to make 
it independent from the reconstructed distance [6T] . 
This method gives good results for multiwavelength 
hologram multiplexing, but may however, increase 
the computational load. Another way to adjust the 
magnification is to reconstruct the hologram with a 
two step algorithm [62]. Each step consists of a 1-FFT 
reconstruction. Let z be the reconstruction distance, 
the two steps (reconstruction at distances zi and Z2) 
are chosen such that z = zi -\- Z2. Here, the mag- 
nification is controlled by the choice of the interme- 
diate reconstruction distance zi. An optimization of 
this approach allows the authors of Ref. [63J to bet- 
ter match physical diffraction, thus obtaining high- 
fidelity reconstruction of magnified holograms. Con- 
trol of reconstruction magnification, shift, and aber- 
ration compensation has also been proposed and real- 
ized using a digital lens, with adjustable parameters, 
in the reconstruction process [64j. 

In the following subsection, we will focus on two al- 
gorithms allowing the adjustment of magnification in 
the reconstruction process and that are based either 
on the convolution [65] or the 1-FFT implemen- 
tation of the Fresnel transform. 

3.C.I. Digital quadratic lens method 

This method is based on the convolution ap- 
proach [591 160]. Prior to reconstruction, the hologram 
is padded to the desired horizon and then multiplied 
by a digital spherical wavefront, acting as a quadratic 
lens, which is defined by: 



C{x) = exp ( 



(15) 



where Rc denotes the curvature radius of C. This cur- 
vature radius can be defined in terms of system mag- 
nification such that: 



jz 



(16) 



Here, 7 is the ratio between the CCD horizon (of 
the padded hologram) and the object physical ex- 
tent. Working with a spherical reconstruction wave- 
front modifies the physical reconstruction distance z 
to z' = jz. Thus, hologram reconstruction can be 
realized by computing the following relation: 



exp {'^z') 
T-'[T{E{x)C{x)}T{h,^{x)}], (17) 



^rec (0 



when working within a "three-FFT" scheme, or 

(18) 

when angular spectrum propagation is considered. 
An alternative method, based on this formalism as- 
sociated with a spatial filtering of the 1-FFT re- 
constructed hologram, allows the reconstruction of 
the local object field with an adjustable magnifica- 
tion [65j. This method makes it possible to limit the 
effect of the reference beam distortions. However, one 
more FFT (two, when the angular spectrum imple- 
mentation is considered) is needed to deal with the 
filtering step. 

3.C.2. Fresnel-Bluestein transform 

This approach is based on a "clever" expansion of 
Eq. (|9|) [66j. In the kernel of the Fourier transform, 
the product 2np is rewritten as 2np = + — 
{p — nf' [67j, such that the discrete Fresnel transform 
can be expressed by: 



^ ^ exp i^^z 

Erec {P) = . , ^ exp 



iXz 
E (n) exp 



-—A^Ax-AOp^ 



N 



n=0 



—Ax {Ax - A^) 

AZ 



X exp 



Xz 



AxA^ {p — n) 



■ (19) 



Let 7 = A^/Ax be the magnification of the recon- 
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Fig. 3: Reconstruction of an hologram for 7 = 
0.8, 1, 2.5 X 7o. (a,c,e) Quadratic lens method. (b,d,f) 
Fresnel-Bluestein method. 

struction. Therefore, Eq. ([T9|) can be rewritten as: 



^rec (P) = A ^ exp 



-^7(1-7) AxV 



N 



exp 



n=0 



-^(l-7)nW 



X exp 



7 (p — n) Ax^ 



(20) 



It should be noted that the magnification 7 is inde- 
pendent of the hologram recording parameters, and 
can be adjusted at will. With this formulation, Eq. 
(f2Q|) is the spatial convolution product of two func- 
tions / and g defined by: 



f{n) = E{n)exp i— {1 - j) n'^ Ax'^ 



(21) 



g{n) = exp (^i^jn'^Ax'^^ . (22) 



The Fresnel-Bluestein reconstruction algorithm can, 
therefore, be summarized as: 



^rec (0 



exp (^z) 



iXz 



exp 



-—7(1-7) AxV 



:F-'[:F{f{x)}:F{g{x)}]. (23) 



Adjustable magnification rendering with the 
quadratic lens method and the Fresnel-Bluestein al- 
gorithm yield to the same results. Let 70 be the 
intrinsic magnification of the 1-FFT reconstruction 
scheme. Hologram rendering of a USAF resolution 
target sector, with 228 line pairs. mm~^ spatial fre- 
quency at element (7-6), with both methods at mag- 
nification 7 = 0.8, 1, 2.5 X 70 is reported in Fig. (|3]). 
Despite these two methods are based on different 
formalisms (single-FFT formalism, and convolution 
based approach), it is here made obvious that both 
methods lead to the same results. 

3.C.3. Aliases and replicas 

Adjustable magnification algorithms make it possible 
either to zoom over details in reconstructed images or 
to reconstruct objects whose dimensions are greater 
than that of the recording device. However, cares are 
to be taken. As a matter of fact, working with recon- 
struction horizons smaller than the object physical 
extend may cause aliases in the reconstruction plane, 
whereas replicas may appear in the opposite situa- 
tion. In other words, considering 70 = Xz/ [NAx'^) 
the intrinsic magnification of the 1-FFT based Fres- 
nel transform implementation. 



7 < 70 



(24) 



will lead to replicas in the reconstructed image, 
whereas choosing 

7 > 70 (25) 

will generate aliases. 

This aspect is illustrated by Fig.|3](a) and|3](c). A 
hologram of the object was reconstructed using the 
adjustable magnification algorithm proposed by Re- 
strepo ^66j (results would have been the same if the 
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Fig. 4: Illustration of alias and replica phenomena, (a) Reconstruction with 7 = 0.5 x 70. (c) Reconstruc- 
tion with 7 = 4 X 70. (d) Same as (a) with replica removal, (f) Same as (c) with aliases filtering. (b,e) 
Reconstruction with 7 = 70. 
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Fig. 5: (a) Synoptics of the anti-alias procedure, (b) 
Replica removal scheme. 



quadratic lens algorithm was considered) . The recon- 
struction with the 1-FFT algorithm (7 = 70) is pro- 



posed on Fig. (jll (b) and (e). It should be noted that 
replicas can be seen in Fig. |4] (a) (here 7 < 70) and 
aliases appear in Fig. [H (c) (for 7 > 70). These un- 
wanted effects degrade the reconstruction quality and 
must be avoided. 

To limit the aliasing effect, when 7 > 70, Hen- 
nelly proposed a filtering scheme which is presented 
in Fig. ([5]) (a) ^68j. Prior to being reconstructed, the 
hologram is multiplied by a chirp function defined as: 



C (x) = exp 



(26) 



The resulting chirp-multiplied hologram is then low- 
pass filtered and finally multiplied by C*, which is the 
complex conjugate of C. The size of the filtering win- 
dow is chosen so as to match the physical extent of 
the reconstruction horizon at distance z. It is there- 
fore possible to reconstruct the hologram using an ad- 
justable magnification algorithm. Benefits of this fil- 
tering approach is illustrated by Figs.lH (c) andlH (f). 
Aliases artifacts are completely removed, thus giving 
a high contrast image of the reconstructed objects. 
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When 7 < 7o, replicas can be seen on the recon- 
structed image. Their removal can be performed by 
the procedure illustrated in Fig. ([5]) (b). The recon- 
structed hologram is cropped in order to keep the 
(70/7) A/" pixels associated with the original object. 
This selection is then zero-padded to the original size 
of the hologram. As can be seen from Figs. 2] (a) and 
m (d) , replicas are completely removed. 

3.D. Fresnelet decomposition 

Fresnelet decomposition was initially proposed by 
Liebling for the reconstruction and processing of dig- 
ital holograms ^69j . This multiresolution scheme finds 
application in a wide variety of domains such as data 
compression [TQ l ITT ] , non linear filtering [72j, wave- 
front retrieving [73] and can be considered in autofo- 
cusing procedures [74]. Fresnelet reconstruction of a 
hologram consists of its decomposition on a basis of 
Fresnel-transformed wavelets. 

Liebling proposed the use of B-splines, which can 
be defined as [69j: 



(x) = * . . . * (x) , 



(27) 



n+l 



where is given by: 



1, 0<x <1 
p^{x) = ^ 1, x = Oorx=l , 
0, otherwise 



(28) 



and the * symbol denotes the convolution product. As 
shown by Unser, B-spline fulfills all the mathematical 
requirements to be used for multiresolution analysis 
of L2 (M) [75J, especially the two-scale relation: 



(29) 



Here h{k) = ih: i ^ | is the binomial filter. 

V ^ / 

B-spline can be used to generate a semi-orthogonal 
wavelet function basis of L2 (M) denoted ip^j^ and de- 
fined as: 



j,kez 



(30) 



where, 



r{^)=Y.9ik)P''{x-k). (31) 
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Fig. 6: Experimental procedure for the holographic 
reconstruction benchmarking 

The filter g (k) is the quadrature mirror filter of h {k). 
Fresnelets basis can be calculated simply by taking 
the Fresnel transform of the B-spline basis. Fresnelet 
bases are therefore defined by: 



(32) 



with 



and 



r{^)=^g{k)^"{x-k), (33) 



/3"(f)=E^(^)^" (^-fc)- 



(34) 



Here 7 is associated with the Fresnel transform. It 
should be noted that computation method chosen for 
the Fresnel transform will affect the Fresnelet trans- 
form results. As a matter of fact, Fresnelet transform 
properties will be the same as those of the chosen 
Fresnel computation scheme {e.g. adjustable, or uni- 
tary magnification). 

4. Application 

In this part, reconstruction of experimental holo- 
grams is performed according to the methods re- 
ported in Section (|3|). Holograms are recorded ac- 
cording to the experimental set-up of Fig. ([6|). Here, 
off- axis interference between the reference and object 
beam are recorded on a 2048 x 2048 pixels CCD sen- 
sor with Ax = 7.4 /im pixel pitch. The object consists 
of an inverted USAF target illuminated with a green 
laser (A = 532 nm). The USAF target is positioned at 
three different distances z^, from the sensor, chosen 
such that < Ax, A^ = Ax, and A^ > Ax, where 
A^ and Ax denote the size of the reconstruction and 
of the CCD sensor respectively. Experimental recon- 
struction are presented hereafter. 
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4. A. Classical reconstruction methods 
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In this section, holograms recorded at Zi such that 
< Ax, A^ = Ax, and A^ > Ax, are recon- 
structed using 1-FFT, angular spectrum propagation, 
and 3-FFT methods. Images of the reconstructed ob- 
jects are proposed in Fig. ([7]). This figure consists of 
a two entry table. In each row, the hologram recon- 
structed when A^ < Ax, A^ = Ax, and A^ > Ax 
are depicted. Each column is associated with the cho- 
sen reconstruction method: 1-FFT, angular spectrum 
propagation, and 3-FFT. It is noticeable that in most 
cases reconstruction result depend on the method 
chosen. In the following section, results obtained are 
detailed row by row. 

4.A.I. Reconstruction for < Ax 
1. 1-FFT 

As seen on Fig. ([7]), aliases are present in the 
reconstructed image of the object. This is due to 
the fact that the 1-FFT implementation of the 



Fresnel transform results in a magnified image 
of the original object. Thus, the reconstructed 
object extend over the limits of the CCD sensor. 



2. Angular spectrum 

As far as A^ < Ax, the reconstructed object 
is well embedded within the CCD sensor hori- 
zon. Angular spectrum method is therefore well 
suited for reconstruction of holograms recorded 
near the CCD sensor. 



3. 3-FFT 

The reconstructed image of the object is embed- 
ded within the CCD sensor. However, replicas 
can be noticed from the reconstructed hologram. 
This is due to the fact that, when the reconstruc- 
tion distance z < NAx/X, the impulse response 
hz is ill-sampled: in this situation, the sampling 
theorem is not verified [76] . 
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4. A. 2. Reconstruction for = Ax 

It can be noted that, for = Ax, the three re- 
construction methods considered give the same re- 
sults. As a matter of fact, in this situation, the re- 
constructed horizon perfectly matches the sensor ar- 
ray extend. In other words, intrinsic magnification of 
1-FFT Fresnel implementation is the same as that of 
convolution approaches. 

4. A. 3. Reconstruction for A^ > Ax 

1. 1-FFT 

As far as the object extent is bigger than the sen- 
sor array dimensions, the 1-FFT implementation 
of the Fresnel transform is appropriate for holo- 
gram reconstruction. 

2. Angular spectrum and 3-FFT 

The fact that these approaches exhibit unitary 
magnification is limiting when dealing with an 
object located far from the sensor. As a mat- 
ter of fact, it can be realized from Fig. (O that 
aliases occur in the reconstructed image of the 
hologram. 

Classical reconstruction methods have been ap- 
plied to experimental holograms recorded at various 
distances from the sensor array. It can be noted that 
each method is valid only within a limited range of 
distances. In the next section, we will give a few word 
about Fresnelet decomposition, and show that its 
reconstruction properties can be modified to match 
each reconstruction method. 

4.B. Fresnelets 

As presented in Section IXD) Fresnelet decomposition 
is similar to a multiscale-wavelet decomposition on 
a Fresnel-transformed base. One appealing feature of 
this decomposition is that the result of the Fresnelet 
reconstruction depends on the method chosen to com- 
pute the Fresnel transform of the wavelet base. 

To illustrate this point, Fresnelet reconstruction of 
the hologram recorded for A^ > Ax is performed. In 
this situation, the 1-FFT method gave good results, 
whereas the 3-FFT reconstruction produces alises. 
Here, the fresnelet bases are calculated with the 1- 
FFT and the 3-FFT method according to Eqs. §^ 




(c) (d) 



Fig. 8: Fresnelet decomposition of the hologram 
recorded for A^ > Ax. (a) Fresnelet coefficients com- 
puted within the 1-FFT scheme, (b) Fresnelet coeffi- 
cients computed within the 3-FFT scheme, (c) Holo- 
gram reconstruction from (a), (d) Hologram recon- 
struction from (b). 

and (|34|) . Decompositions of the test hologram on the 
two calculated fresnelet bases are proposed in Fig. ([8]) 
(a) and (b) respectively. It should be noticed that the 
computation scheme chosen strongly affects the cal- 
culated coefficients. Therefore, properties of the fres- 
nelet decomposition reconstruction depends on the 
method chosen to calculate the fresnelet base func- 
tions. This aspect is pointed out by Fig. (|8]) (c) and 
(d). Here, hologram reconstruction from the Fresnelet 
coefficients depicted in Fig. (jSj) (a) and (b) is real- 
ized. These reconstructions are similar to the one ob- 
tained with classical methods (See A(f > Ax in Fig. 
([7j) for comparison). Thus, for A^ > Ax, the single- 
FFT method will be more reliable than the 3-FFT 
scheme. 

5. Conclusion 

We have proposed an overview of holographic recon- 
struction methods with application to off-axis inten- 
sity hologram treatment. Intrinsic properties and lim- 
itations of classical methods have been investigated 
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and applicability ranges have been stated with the 
reconstruction of experimental holograms. It can be 
noted that the choice of the reconstruction method 
will be driven by the hologram recording conditions. 
For instance, when dealing with far and extended ob- 
jects, 1-FFT algorithm will be the most appropri- 
ate, whereas convolution approaches will be suited 
for the reconstruction of small objects located near 
the sensor array. Adjustable magnification methods 
have been presented and can be viewed as a way 
to overcome limitations of the classical reconstruc- 
tion schemes and allow the reconstruction result to 
be independent from the chosen scheme. Neverthe- 
less, cares are to be taken in order to limit aliases and 
replicas when working with high or low magnification. 
Finally Fresnelet reconstruction of the holograms has 
been performed using the fact that the Fresnelet de- 
composition base depends on the method chosen to 
compute the Fresnel transform. This method can be 
considered for filtering or image compression when 
computational load is not a critical issue. 

References 

1. D. Gabor, "A new microscopic principle," Nature 
161, 777-778 (1948). 

2. D. Gabor, "Microscopy by reconstructed wave- 
fronts," Proc. Roy. Soc. A. 197, 454-487 (1949). 

3. E. N. Leith, J. Upatnieks, "Reconstructed wave- 
fronts and communication theory," J. Opt. Soc. 
Am. 52, 1123-1130 (1962). 

4. J.W. Goodman, R.W. Lawrence, "Digital im- 
age formation from electronically detected holo- 
grams," Appl. Phys. Lett. 11, 77-79 (1967). 

5. M.A. Kronrod, N.S. Merzlyakov, L.P. 
Yaroslavsky, "Reconstruction of holograms 
with a computer," Sov. Phys. -Tech. Phys. 17, 
419-420 (1972). 

6. L.P. Yaroslvsky, N.S. Merzlyakov, Methods of 
digital holography. Consultant Bureau, New 
York, (1980). 

7. L. Onural, P. D. Scott, "Digital recording of in- 
line holograms," Opt. Eng. 26, 1124-1132 (1987). 

8. G. Liu, P. D. Scott, "Phase retrieval and twin- 
image elimination for in-line Fresnel holograms," 
J. Opt. Soc. Am. A 4, 159-165 (1987). 



9. J.R. Fienup, "Phase retrieval algorithms, a com- 
parison," Appl. Opt. 21, 2758-2769 (1982). 

10. L. Onural, M. T. Ozgen, "Extraction of three- 
dimensional object-location information directly 
from in-line holograms using Wigner analysis," 
J. Opt. Soc. Am. A 9, 252-260 (1992). 

11. U. Schnars, W. Jiiptner, "Direct recording of 
holograms by a CCD target and numerical re- 
construction," Appl. Opt. 33, 179-181 (1994). 

12. A. Lozano, J. Kostas, J. Soria, "Use of holog- 
raphy in particle image velocimetry measure- 
ments of a swirling flow," Exp. Fluids 27, 251- 
261 (1999). 

13. H. Meng, G. Pan, Y. Pu, S.H. Woodward, "Holo- 
graphic particle image velocimetry: from film to 
digital recording," Meas. Sci. Technol. 15, 673- 
685 (2004). 

14. Y. Pu, H. Meng, "Four-dimensional dynamic 
flow measurement by holographic particle image 
velocimetry," Appl. Opt. 44, 7697-7708 (2005). 

15. F. Dubois, N. Callens, C. Yourassowsky, M. 
Hoyos, P. Kurowski, O. Monnom, "Digital holo- 
graphic microscopy with reduced spatial coher- 
ence for three-dimensional particle flow analy- 
sis," Appl. Opt. 45, 864-871 (2006). 

16. M. Atlan, M. Gross, "Laser Doppler imaging, 
revisited," Rev. Sci. Instruments 77, 116103 
(2006). 

17. J-M. Desse, P. Picart, P. Tankam, "Digital three- 
color holographic interferometry for flow analy- 
sis," Opt. Express 16, 5471-5480 (2008). 

18. N. Verrier, S. Coetmellec, M. Brunei, D. Le- 
brun, "Digital in-line holography in thick optical 
systems: application to visualization in pipes," 
Appl. Opt. 47, 4147-4157 (2008). 

19. N. Verrier, S. Coetmellec, M. Brunei, D. Lebrun, 
"Determination of 3D-region of interest using 
digital in-line holography with astigmatic Gaus- 
sian beams," J. Europ. Opt. Soc. Rap. Public. 4, 
09038 (2009). 

20. N. Verrier, C. Remacha, M. Brunei, D. Lebrun, 
S. Coetmellec, "Micropipe flow visualization us- 
ing digital in-line holographic microscopy," Opt. 
Express 18, 7807-7819 (2010). 



11 



21. S. Schedin, G. Pedrini, H.J. Tiziani, "Pulsed dig- 
ital holography for deformation raeasurements 
on biological tissues," Appl. Opt. 39, 2853-2857 
(2000). 

22. M.K. Kim, "Tomographic three-dimensional 
imaging of a biological specimen using 
wavelength-scanning digital interference holog- 
raphy," Opt. Express 7, 305-310 (2000). 

23. F. Charriere, N. Pavilion, T. Colomb, C. De- 
peursinge, T.J. Heger, E.A.D. Mitchell, P. Mar- 
quet, B. Rappaz, "Living specimen tomography 
by digital holographic microscopy: morphometry 
of testate amoeba," Opt. Express 14, 7005-7013 
(2006). 

24. W. Xu, M. H. Jericho, I. A. Meinertzhagen, H. J. 
Kreuzer, "Digital in-line holography for biologi- 
cal applications," PNAS 98, 11301-11305 (2001). 

25. B. Kemper, G. von Bally, "Digital holographic 
microscopy for live cell applications and technical 
inspection," Appl. Opt. 47, A52-A61 (1994). 

26. M. Simonutti, M. Paques, J. A. Sahel, M. Gross, 
B. Samson, C. Magnain, and M. Atlan, "Holo- 
graphic laser Doppler ophthalmoscopy," Opt. 
Lett. 35, 1941-1943 (2010). 

27. R.L. Powell, K.A. Stetson, "Interferometric vi- 
bration analysis by wavefront reconstruction," J. 
Opt. Soc. Am. 55, 1593-1597 (1965). 

28. C.C. Aleksoff, "Temporally modulated hologra- 
phy," Appl. Opt. 10, 1329-1341 (1971). 

29. F. Zhang, J.D.R. Valera, I. Yamaguchi, M. 
Yokota, G. Mills, "Vibration analysis by phase 
shifting digital holography," Opt. Rev. 11, 297- 
299 (2004). 

30. U. lemma, L. Morino, M. Diez, "Digital holog- 
raphy and Karhunen-Loeve decomposition for 
the modal analysis of two-dimensional vibrat- 
ing structures," J. Sound Vibration 291, 107-131 
(2006). 

31. P. Picart, J. Leval, D. Mounier, S. Gougeon, 
"Some opportunities for vibration analysis with 
time averaging in digital Fresnel holography," 
Appl. Opt. 44, 337-343 (2005). 

32. J. Leval, P. Picart, J-P. Boileau, J-C. Pascal, 
" Full- field vibrometry with digital Fresnel holog- 
raphy," Appl. Opt. 44, 5763-5772 (2005). 



33. D. Borza, "Mechanical vibration measurement 
by high-resolution time-averaged digital hologra- 
phy," Meas. Sci. Technol. 16, 1853-1864 (2005). 

34. A. Asundi, V.R. Singh, "Time- aver aged in-line 
digital holographic interferometry for vibration 
analysis," Appl. Opt. 45, 2391-2395 (2006). 

35. F. Joud, F. Lanoe, M. Atlan, J. Hare, M. Gross, 
"Imaging a vibrating object by sideband digital 
holography," Opt. Express 17, 2774-2779 (2009). 

36. I. Yamaguchi, K. Yamamoto, G.A. Mills, M. 
Yokota, "Image reconstruction only by phase 
data in phase-shifting holography," Appl. Opt. 
45, 975-983 (2006). 

37. K. Matsushima, "Shifted angular spectrum 
method for off- axis numerical propagation," Opt. 
Express 18, 18453-18463 (2010). 

38. D. Lebrun, A. Benkouider, S. Coetmellec, M. 
Malek, "Particle field digital holographic recon- 
struction in arbitrary tilted planes," Opt. Ex- 
press 11, 224-229 (2003). 

39. S. De Nicola, A. Finizio, G. Pierattini, P. Fer- 
raro, D. Alfieri, "Angular spectrum method with 
correction of anamorphism for numerical recon- 
struction on tilted planes," Opt. Express 13, 
9935-9940 (2005). 

40. N. Verrier, S. Coetmellec, M. Brunei, D. Le- 
brun, A.J.E.M. Janssen, "Digital in-line holog- 
raphy with an elliptical, astigmatic, Gaussian 
beam: wide angle reconstruction," J. Opt. Soc. 
Am. A 25, 1459-1466 (2008). 

41. F. Soulez, L. Denis, C. Fournier, E. Thiebaut, 

C. Goepfert, "Inverse-problem approach for par- 
ticle digital holography: accurate location based 
on local optimization," J. Opt. Soc. Am. A 24, 
1164-1171 (2007). 

42. F. Soulez, L. Denis, E. Thiebaut, C. Fournier, C. 
Goepfert, "Inverse-problem approach in particle 
digital holography: out-of-field particle detection 
made possible," J. Opt. Soc. Am. A 24, 3708- 
3716 (2007). 

43. L. Denis, D. Lorenz, E. Thiebaut, C. Fournier, 

D. Trede, "Inline hologram reconstruction with 
sparsity constraints," Opt. Lett 34, 3475-3477 
(2009). 

44. M. Marim, E. Angelini, J-C. Olivo-Marin, M. 



12 



Atlan, "Off-axis compressed holographic rai- 
croscopy in low-hght conditions," Opt. Lett 36, 
79-81 (2011). 

45. T. Shimobaba, Y. Sato, J. Miura, M. Take- 
nouchi, T. Ito, "Real-time digital holographic mi- 
croscopy using the graphic processing unit," Opt. 
Express 16, 11776-11781 (2008). 

46. L. Ahrenberg, A.J. Page, B.M. Hennelly, J.B. 
McDonald, T.J. Naughton, "Using commodity 
graphics hardware for real-time digital hologram 
view-reconstruction," IEEE J. Display Technol. 
5, 111-119 (2009). 

47. T. Shimobaba, N. Masuda, Y. Ichihashi, T. Ito, 
"Real-time digital holographic microscopy ob- 
servable in multi-view and multi-resolution," J. 
Opt. 12, 065402 (2010). 

48. B. Samson, F. Verpillat, M. Gross, M. Atlan, 
"Video-rate wide-field laser vibrometry by het- 
erodyne holography," Opt. Lett. 36, 1449-1451 
(2011). 

49. I. Yamaguchi, T. Zhang, "Phase-shifting digital 
holography," Opt. Lett. 22, 1268-1270 (1997). 

50. M. Atlan, M. Gross, E. Absil, "Accurate phase- 
shifting digital interferometry, " Opt. Lett. 32, 
1456-1458 (2007). 

51. J.W. Goodman, Introduction to Fourier Optics, 
Roberts and Company, third edition, (2005). 

52. E. N. Leith, J. Upatnieks, "Wave front recon- 
struction with continuous-tone objects," J. Opt. 
Soc. Am. 53, 1377-1381 (1963). 

53. E. Cuche, P. Marquet, C. Depeursinge, "Spatial 
filtering for zero-order and twin-image elimina- 
tion in digital off-axis holography," Appl. Opt. 

39, 4070-4075 (2000). 

54. I. Yamaguchi, J-i. Kato, S. Ohta, J. Mizuno, "Im- 
age formation in phase-shifting digital hologra- 
phy and applications to microscopy," Appl. Opt. 

40, 6177-6186 (2001). 

55. M. Born, E. Wolf, Principles of optics, Cam- 
bridge University Press, seventh edition, (1999). 

56. J. W. Cooley and J. W. Tukey, "An Algorithm 
for the machine computation of complex Fourier 
series," Math. comp. 19, 297-301 (1965). 

57. U. Schnars, W. Jiiptner, "Digital recording and 



numerical reconstruction of holograms," Meas. 
Sci. Technol. 13, R85-R101 (2002). 

58. L. Yu, M. K. Kim, "Wavelength-scanning digital 
interference holography for tomographic three- 
dimensional imaging by use of the angular spec- 
trum method," Opt. Lett. 30, 2092-2094 (2005). 

59. J-c. Li, P. Tankam, Z-j. Peng, P. Picart, "Digital 
holographic reconstruction of large object using 
a convolution approach and adjustable magnifi- 
cation," Opt. Lett. 34, 572-574 (2009). 

60. P. Picart, P. Tankam, D. Mounier, Z-j. Peng, J- 
c. Li, "Spatial bandwidth extended reconstruc- 
tion for digital color Fresnel holograms," Opt. 
Express 17, 9145-9156 (2009). 

61. P. Ferraro, S. De Nicola, G. Coppola, A. Finizio, 
D. Alfieri, G. Pierattini, "Controlling image size 
as a function of distance and wavelength in 
Fresnel-transform reconstruction of digital holo- 
grams," Opt. Lett. 26, 854-856 (2004). 

62. F. Zhang, I. Yamaguchi, L.P. Yaroslavsky, "Al- 
gorithm for reconstruction of digital holograms 
with adjustable magnification," Opt. Lett. 29, 
1668-1670 (2004). 

63. D. Wang, J. Zhao, F. Zhang, G. Pedrini, W. 
Osten, "High-fidelity numerical realization of 
multiple-step Fresnel propagation for the recon- 
struction of digital holograms," Appl. Opt. 47, 
D12-D20 (2008). 

64. T. Colomb, F. Montfort, J. Kiihn, N. Aspert, E. 
Cuche, A. Marian, F. Charriere, S. Bourquin, P. 
Marquet, C. Depeursinge, "Numerical paramet- 
ric lens for shifting, magnification, and complete 
aberration compensation in digital holographic 
microscopy," J. Opt. Soc. Am. A. 23, 3177-3190 
(2006). 

65. J-c. Li, Z-j. Peng, P. Tankam, Q-h. Song, P. Pi- 
cart, "Digital holographic reconstruction of a lo- 
cal object field using an adjustable magnifica- 
tion," J. Opt. Soc. Am. A 28, 1291-1296 (2011). 

66. J. F. Restrepo, J. Garcia-Sucerquia, "Magnified 
reconstruction of digitally recorded holograms 
by Fresnel-Bluestein transform," Appl. Opt. 49, 
6430-6435 (2010). 

67. L. Bleustein, "Linear filtering approach to the 
computation of the discrete Fourier transform," 



13 



IEEE Trans. Audio Electroacoust. 18, 451-455 
(1970). 

68. B. Hennelly, D. Kelly, N. Pandey, D. Monaghan, 
"Zooming algorithms for digital holography," J. 
Phys.: Conf. Series 206, 012027 (2010). 

69. M. Liebling, T. Blu, M. Unser, "Fresnelet: new 
multiresolution wavelet bases for digital holog- 
raphy," IEEE Trans. Image Process. 12, 29-43 
(2003). 

70. E. Darakis, J. J. Soraghan, "Use of Fresnelets for 
phase-shift digital hologram compression," IEEE 
Trans. Image Process. 15, 3804-3811 (2006). 

71. E. Darakis, T. J. Naughton, J. J. Soraghan, 
"Compression defect in different reconstructions 
from phase-shifting digital holographic data," 
Appl. Opt. 46, 4579-4586 (2007). 

72. M. Liebling, T. Blu, M. Unser, "Non-Linear Fres- 
nelet Approximation for Interference Term Sup- 
pression in Digital Holography", Proceedings of 



the SPIE Conference on Mathematical Imaging: 
Wavelet Applications in Signal and Image Pro- 
cessing X, San Diego CA, USA, 5207, 553-559 

(2003) . 

73. M. Liebling, T. Blu, M. Unser, " Complex- wave 
retrieval from a single off-axis hologram," J. Opt. 
Soc. Am. A 21, 367-377 (2004). 

74. M. Liebling, M. Unser, "Autofocus for digital 
Fresnel holograms by use of a Fresnelet-sparsity 
criterion," J. Opt. Soc. Am. A 21, 2424-2430 

(2004) . 

75. M. Unser, A. Aldroubi, M. Eden, "A family of 
polynomial spline wavelet transforms," Signal 
Process. 30, 141-162 (1993). 

76. L. Onural, "Sampling of the diffraction field," 
Appl. Opt. 15, 3804-3811 (2007). 



14 



